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Abstract — In  multicarrier  systems,  when  the  order  of  a  channel 
impulse  response  is  larger  than  the  length  of  the  cyclic  prefix  (CP), 
there  is  a  significant  performance  degradation  due  to  interblock 
interference  (IBI).  This  paper  proposes  a  blind-channel  short¬ 
ening  method  in  which  the  equalizer  parameter  vector  is  formed 
by  the  noise  subspace  of  the  received  signal  correlation  matrix  so 
that  the  output  power  is  maximized.  The  proposed  method  can 
not  only  shorten  the  effective  channel  impulse  response  to  within 
the  CP  length  but  also  maximize  the  output  signal-to-interfer- 
ence-and-noise  ratio  while  eliminating  the  IBI.  We  point  out  that  the 
performance  depends  on  the  choice  of  a  decision  delay  and  propose 
a  simple  method  for  determining  the  appropriate  delay.  We  propose 
both  a  batch  algorithm  and  an  adaptive  algorithm  and  show  by 
simulation  that  they  are  superior  to  the  conventional  algorithms. 

Index  Terms — Blind  equalization,  cyclic  prefix  (CP),  orthogonal 
frequency-division  multiplexing  (OFDM),  time-domain  equalizer. 


I.  Introduction 

MULTICARRIER  modulation,  such  as  orthogonal  fre¬ 
quency-division  multiplexing  (OFDM),  is  an  attractive 
technique  for  high-speed  data  transmission  and  holds  great 
potential  for  audio/video  broadcasting,  wireless  local  area 
networks,  and  future  wideband  cellular  systems.  In  multicarrier 
systems,  a  cyclic  prefix  (CP)  is  inserted  between  data  blocks 
to  prevent  interblock  interference  (IBI)  caused  by  a  dispersive 
channel.  When  the  channel  impulse  response  is  extremely 
long,  a  long  CP  is  not  always  available  since  it  reduces  the 
transmission  efficiency.  If  the  channel  impulse  response  is 
longer  than  the  length  of  CP  plus  one,  IBI  remains.  There 
have  been  some  works  to  cope  with  the  IBI  in  multicarrier 
applications  such  as  Digital  Subscriber  Line  (DSL)  [1],  [2], 
the  Digital  Terrestrial  Television  Broadcasting  (DTTB)  using 
Single  Frequency  Network  (SFN)  [3],  5-GHz  wireless  LAN 
[4],  and  WiMAX  [5],  Channel  shortening  filters  are  known 
to  be  relatively  simple  and  effective  techniques  [1],  [2],  They 
shorten  the  impulse  response  of  the  effective  channel,  which  is 
a  convolution  of  a  channel  and  a  linear  equalizer,  to  within  the 
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CP  length  plus  one.  A  number  of  channel  shortening  methods 
have  been  proposed  in  the  past  decade,  but  most  of  them  require 
training  sequences  that  reduce  channel  throughput. 

The  recently  proposed  blind-channel  shortening  techniques 
[6]— [  12]  are  attractive  because  they  do  not  require  the  transmis¬ 
sion  of  training  sequences.  The  sum-squared  autocorrelation 
minimization  (SAM)  algorithm  [7]  is  an  adaptive  algorithm  that 
tries  to  minimize  the  sum  of  the  squares  of  autocorrelations  of 
all  lags  greater  than  the  desired  window.  The  SAM  cost  function 
has  undesired  local  minima,  though,  and  its  global  convergence 
is  not  assured.  As  demonstrated  in  [10],  it  may  fail  to  shorten 
a  channel  impulse  response.  The  single  lag  autocorrelation 
minimization  (SLAM)  algorithm  [8]  minimizes  only  a  single 
autocorrelation  term.  It  can  be  regarded  as  a  simplified  version 
of  the  SAM  algorithm  and  its  cost  function  also  has  undesired 
local  minima.  An  adaptive  algorithm  exploiting  the  existence 
of  null  subcarriers  [9]  tries  to  force  the  receiver’s  fast  Fourier 
transform  (FFT)  output  corresponding  to  the  null  subcarriers  to 
be  zero.  Although  this  algorithm  requires  frequency-hopping 
transmission  for  global  convergence,  frequency-hopping  null 
tones  are  impractical.  Second-order  statistical  methods  [10] 
search  for  an  equalizer  parameter  vector  in  the  null  space 
of  correlation  matrices,  but  their  computational  complexity 
is  high.  Unlike  the  other  methods,  they  can  be  implemented 
only  by  batch  processing  algorithms.  To  deal  with  both  static 
and  time  varying  environments,  an  approach  should  be  imple¬ 
mented  by  both  batch  and  adaptive  algorithms.  The  multicarrier 
equalization  by  restoration  of  redundancy  (MERRY)  algorithm 
in  [11]  is  a  simple  blind  algorithm,  is  globally  convergent,  and 
can  also  be  implemented  by  a  batch  algorithm. 

This  paper  shows  that  although  the  MERRY  algorithms 
using  over-sampled  channel  outputs  can  achieve  perfect 
channel  shortening,  which  means  that  IBI  can  be  canceled 
completely,  their  output  signal-to-interference-and-noise  ratio 
(SINR)  can  be  unsatisfactory.  We  therefore  propose  a  method 
that  maximizes  the  output  SINR  while  eliminating  the  IBI  and 
that  can  be  implemented  by  both  a  batch  algorighm  and  an 
adaptive  algorithm.  We  also  consider  the  choice  of  a  decision 
delay  that  affects  the  SINR  performance. 

This  paper  is  organized  as  follows.  In  Section  II,  we  first 
describe  a  communication  model  using  a  channel  shortening 
equalizer  and  then  describe  the  MERRY  algorithms.  Section  III 
describes  our  method  for  maximizing  the  output  SINR  while 
eliminating  the  IBI  and  also  presents  a  simple  method  for 
choosing  an  appropriate  decision  delay  .  Section  IV  presents 
batch  and  adaptive  algorithms  implementing  our  methods. 
Section  V  shows  simulation  results,  and  Section  VI  concludes 
the  paper  by  summarizing  it  briefly. 
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Fig.  1.  Baseband  multicarrier  system  model. 
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II.  Problem  Formulation 
A.  Communication  Model 

We  consider  a  baseband  multicarrier  communication  system 
using  the  CP  (Fig.  1).  The  M-ary-QAM  data  sequences  { ,s /.. } 
generated  at  a  rate  1/TS  are  i.i.d.  with  variance  a; .  The  nth  data 
block  of  length  N  is 

Sn  =  [finiV  SnN+1  "•  S(n+l)lV-l]  0) 

where  the  superscript  T  represents  the  transpose  of  a  matrix. 
Applying  an  inverse  fast  Fourier  transform  (IFFT)  to  s„  pro¬ 
duces  a  time-domain  block  vector  of  length  N 

x„  =  [x0[n]  xi[n\  ■  ■  ■  xN-i[n]]T  =  Fs„  (2) 

where  F  is  an  N -dimensional  IFFT  matrix  that  has  a  unitary 
property  FFff  =  I/v,  where  the  superscript  H  represents  the 
conjugate  transpose  of  a  matrix  and  Ijy  denotes  an  identity  ma¬ 
trix  of  size  N.  A  CP  of  length  P  is  added  to  the  beginning  of 
the  block  x„  to  form  the  nth  multicarrier  block  u„  of  length 
Q  =  N  +  P 


^ nQ 

-xN-P\n\ 

A_ 

UnQ+P- 1 

®jy_i[n] 

U71Q+P 

x0  [n] 

_"{n  r  1  )<«? —  1  _ 

Ji 

■  7 

Obviously,  the  transmitted  signal  uj.  is  cyclostationary.  The  re¬ 
ceived  signal  is  sampled  at  t  =  kTs /p,  where  the  integer  p  is  the 
over-sampling  factor.  The  channel  output  samples  can  be  then 
expressed  as  [10] 


—  ^  ^  unhk—np  T  wk 


where  is  the  sampled  impulse  response  of  the  original 
channel,  which  includes  transmitter  and  receiver  filters  as  well 
as  the  physical  channel,  and  the  AWGN  Wk  is  stationary  and 
independent  of  data  sequences  ,s,;.  Assume  that  h(t)  has  joint 
finite  support  [0,  (M  +  1)TS),  where  the  integer  M  + 1  is  length 


of  the  channel  impulse  response.  Thus,  the  channel  output 
samples  are 


=  ^2  h [i\uk-i  +  wfc 


-r(k+i)P-i  ■ 


h[i]  —  \jHp  hip+ 1  *  *  *  ^(i+lip— l]  (6) 

W k  —  [tFfcp  I  *  *  '  /|fF(/,'+l  )p—  \  ]  (7) 

The  input  to  an  equalizer  consists  of  channel  outputs  sampled 
over  an  L  symbol  period.  Then,  let  Lp  be  the  number  of  sam¬ 
pled  channel  outputs  to  be  collected  in  a  block.  We  define  the 
following  notations: 

u[k]  =  [uk  itk-i  ■  ■  ■  i'k  (8) 

=  [wfc  wfc- 1  •  •  •  w  1_l+i]T ■  (9) 

Moreover,  we  form  an  Lp  x  (L  +  M)  block  Toeplitz  matrix 

-h[0]  h[l]  •••  h  [M\  0  •••  0  - 

0  h[0]  h[l]  •••  h  [M] 

0 

.0  •••  0  h[0]  h[l]  •••  h  [M]. 


A  sampled-channel  output  signal  vector  of  length  Lp  can  be 
expressed  as 


r[fc]  = 


=  Hu[fc]  +  w  [k]. 


-Xk-L+l  - 


We  let  the  parameter  vector  of  a  channel  shortening  equalizer 
t*e  g  =  [5o  9i  •  •  •  9Lp-iY ■  The  output  of  the  equalizer,  known 
as  time -domain  equalizer(TEQ),  can  be  described  by 

yk  =  gHr[k ]  =  gH  (Huffc]  +  w[k\) .  (12) 
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Let  the  effective  channel  impulse  response  be  denoted  by  c  = 
[co  Ci  ■  ■  •  Cl+m- i]T  =  H^g.  The  design  of  a  perfect  channel 
shortening  equalizer  can  be  expressed  as  choosing  g  to  force  all 
samples  c,  except  for  preselected  part  of  P+ 1  consecutive  sam¬ 
ples,  to  be  zero.  Thus,  the  perfect  channel  shortening  equalizer 
can  completely  eliminate  IBI.  Our  immediate  goal  is  to  find  the 
vector  g  performing  the  perfect  channel  shortening  from  only 
the  channel  output. 


as  an  uncorrelated  sequence  as  long  as  L  <  N  —  M,  we  can 
express  the  correlation  matrix  Rq  as 

R0  =  2<72HibiHjbi  +  2aflLp  (16) 

and  thus  the  MERRY  cost  function  becomes 

Jo(g)  =  2<7s2||cibi||2  +  2al  (17) 


B.  MERRY  Algorithms 

MERRY  algorithms  [11]  minimize  the  cost  function  given  by 

■h  (g)  =  E  [\ijp-i+s  -  yQ-i+6\2]  ,  s.t.  ||g||  =  1  (13) 

where  6  6  {0, 1, . . . ,  Q  —  1}  is  the  decision  delay  parameter  and 
E\f\  represents  the  ensemble  average.  We  refer  to  this  function 
to  as  the  MERRY  cost  function. 

The  original  MERRY  algorithm  in  [11]  is  a  stochastic  gra¬ 
dient  algorithm  of  the  MERRY  cost  function.  The  adaptive 
MERRY  algorithm  is  given  by 

f  [k]  =r[(k -1)Q  +  P-1  +  S]~  r[(fc  -  l)Q  +  Q  -  1  +  8] 
e[fc]  =gHr[fc] 

=  g[fc  ~  1]  ~  ^g*[fc]f[fc] 

\\g[k~  1]  —  Ater*[fe]f[fe]|| 

where  //,  is  the  step  gain  and  the  superscript  *  represents  complex 
conjugation.  When  p  =  1,  the  adaptive  MERRY  algorithm  is 
global  convergent  [11]. 

The  MERRY  cost  function  can  be  rewritten  as 


Ms)  =  SHMg;  S.t.  ||g||  =  1 

where  R^  is  an  Lp  X  Lp  matrix  defined  by 


(14) 


R  6  =  E 


(r  [P -l  +  g\-T[Q-l  +  «]) 


(r[P-l  +  8]-r[Q-l  +  6]f  .  (15) 


The  parameter  vector  minimizing  the  MERRY  cost  function  is 
given  by  an  eigenvector  corresponding  to  the  minimum  eigen¬ 
value  of  R,. .  In  the  batch  MERRY  algorithm,  the  correlation 
matrix  R*  is  estimated  by  time  averaging,  and  then  the  param¬ 
eter  vector  is  determined  by  the  eigendecomposition  of  the  es¬ 
timated  correlation  matrix. 

Let  us  look  further  into  the  MERRY  cost  function.  We  let  c,]es 
be  the  desired  signal  component  which  consists  of  P  samples 
of  the  impulse  response  of  the  total  system  c  and  let  c;b;  be  the 
IBI  component  that  is  the  remaining  part  of  c.  Let  Hdes  be  a 
matrix  consisting  of  P  columns  of  H  such  that  Cdes  =  H/^sg, 
and  let  H;bi  be  the  remaining  part  of  H  such  that  c;b;  =  H^g- 
When  <5  =  0,  Hdes  is  the  first  P  column  vectors  of  H  and  H;b] 
is  the  remaining  part.  That  is,  H  =  [Hdes  H,bl].  Taking  into 
account  that  the  channel  input  sequence  {rq,}  can  be  regarded 


where  <r,2  =  E  [|wfc|2].  This  implies  that  the  minimization  of 
the  MERRY  cost  function  minimizes  both  the  IBI  and  the  noise 
component. 

The  MERRY  cost  function  is  reasonable,  but  its  performance 
is  limited.  When  the  noise  is  absent,  the  correlation  matrix  Rq 
is  semi-positive  definite,  Rq  >  0,  from  its  definition.  Moreover, 
when  p  =  1,  Ro  has  full  rank  because 

rank(Ro)  =  rank(H;b;)  =  L. 

Thus,  Ro  is  positive  definite.  Then,  since  Jo  (g)  >  0  for  any 
g  0,  1 1 c;b; 1 1  ^  0.  This  implies  that  the  equalizer  with  p  =  1 
using  the  MERRY  algorithms  cannot  shorten  the  channel  per¬ 
fectly  even  in  the  absence  of  the  noise. 

III.  Perfect  Blind-Channel  Shortening 

A.  Shortenability 

Let  us  consider  the  shortenability  of  the  MERRY  approach 
using  over-sampled  channel  outputs.  That  is,  let  us  consider  it 
when  p  >  2.  For  simplicity  assume  for  the  moment  that  <*>  =  0. 
Also  assume  that  H  has  full  column  rank,  which  is  a  common 
assumption  in  blind  equalization  and  identification  [16].  We 
then  get  the  following  result. 

Proposition  1:  When  p  >  2  and  H  has  full  column  rank, 
the  IBI  component  can  be  canceled  completely  by  minimizing 
Jo  (g)- 

Proof:  Let  Ao  >  Ab  >  •  •  •  >  \lp-i  be  the  eigenvalues 
of  Rq  and  g;,  i  =  0, . . . ,  Lp  —  1,  be  the  corresponding  eigen¬ 
vectors.  From  (14),  the  minimization  of  .Jq  (g)  under  the  norm 
constraint  ||g||  =  1  can  be  done  by  using  g  =  gpp_i.  Since 
H;bi  has  full  column  rank,  the  minimum  eigenvalue  of  R(,  is 
Aip_i  =  2<j,2  .  And  since  Jo(g;)  =  A i,  the  minimum  of  J0  is 
Jo(gLp-i)  =  2<t,2  .  From  (17),  this  implies  that  ||cibi||2  =  0. 
That  is,  it  implies  that  the  IBI  is  canceled  completely.  ■ 

Shortenability  issues  for  more  general  cases  have  been  dis¬ 
cussed  in  [13]. 

B.  SINR  Maximization 

Although  both  the  IBI  and  the  noise  power  can  be  minimized 
by  minimizing  the  MERRY  cost  function,  this  does  not  mean 
that  the  desired  signal  component  is  not  reduced.  As  shown 
later  in  simulation  results,  a  drawback  of  the  MERRY  algorithm 
using  over-sampled  channel  outputs  is  that  it  can  weaken  the  de¬ 
sired  signal  component,  i.e.,  noise  enhancement.  Our  simulation 
result  shown  in  Section  V  suggests  that  when  the  signal-to-noise 
ratio  (SNR)  is  low  this  problem  is  more  serious  when  p  >2  than 
when  p  =  1.  To  overcome  this  problem  we  propose  a  method 
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that  maximizes  the  output  SINR  while  completely  removing  the 

mi. 

Since  the  rank  of  Hibi  is  L  +  M  —  P,  there  are  Lp  —  (L  + 
M  —  P)  eigenvectors  g;  spanning  the  noise  subspace  of  Rq. 
The  following  linear  combination  of  the  q  eigenvectors  is  used 
as  an  equalizer  parameter  vector 

g  =  gLp-iao  +  gLp-2ai  H - b  gLP-qaq -1  =  Gqaq  (18) 

where 

Gq  =\gLp-l  ■  ■  ■  gLp-q]  (19) 

aq  =  [«0  •  •  •  aq- 1]  .  (20) 

The  coefficient  vector  a,/  is  chosen  so  that  the  output  power  is 
maximized,  and  the  cost  function  is  given  by 


Then,  since 

SINR  =  J(ag)~CT^ 

the  SINR  can  be  maximized  by  maximizing  J  ( aq ).  ■ 

We  can  expect  the  SINR  to  increase  as  the  number  q  of 
the  combined  vectors  increases.  We  then  obtain  the  following 
proposition. 

Proposition  3:  Denote  by  Sq  the  maximum  output  power  ob¬ 
tained  by  optimally  combining  q  eigenvectors.  Then  Sq  >  S% 
for  0  <  i  <  q  as  long  as  q  <  Lp  —  (L  +  M  —  P). 

Proof:  Let  ago  be  the  coefficient  vector  that  maximizes  the 
output  power.  That  is 

ag0  =  arg  max  J  (ag) 

IM1=1 

and  we  have 


J(aq)  =E[\yk\2]  ,  s.t.  ||aj  =  1 

=  afGfR,.Ggag  (21) 

where  R r  =  E  [r[/i,-]rH[&]] .  When  q  <  Lp  —  (L  +  M  —  P), 
H^jGgag  =  0  for  any  ag,  and  the  cost  function  can  thus  be 
rewritten  as 

J(aq)  =  us2afGfHdesHfesGgag  +  cr  2 . 

The  first  term  corresponds  to  the  desired  signal  power  and 
the  second  term  corresponds  to  the  noise  power.  Maximizing 
J  (aq)  enables  the  desired  signal  power  to  be  maximized 
without  changing  the  noise  power.  The  basic  idea  of  the 
proposed  method  is  that  IBI  is  canceled  by  choosing  g  from 
the  noise  subspace  of  Rq  and  the  output  SINR  is  maximized 
by  maximizing  the  output  power.  We  obtain  the  following 
proposition. 

Proposition  2:  The  equalizer  parameter  vector  is  given  by 
(18).  Let  u[k]  =  [life  uk-i  ■  ■  ■  Uk-P+i]T ■  When  q  <  Lp  — 
(L  +  M  —  P),  the  output  SINR  defined  by 


SINR  = 


E  [lcLuMI2! 

E  [I \yk  ~  cfesu[fc]|2' 


(22) 


is  maximized  by  choosing  the  coefficient  vector  aq  so  that  the 
output  power  J  (aq)  is  maximized. 

Proof:  Since  {up }  and  { wp. )  are  uncorrelated  sequences 
and  are  uncorrelated  with  each  other  and  ||g||  =  1,  the  SINR 
can  be  rewritten  as 


SINR  = 


^IICdesH2 

o-f||cibi||2  +  oir 


Since  g  is  chosen  from  the  noise  subspace  of  Ro,  IBI  is  canceled 
completely  i.e.,  ||cibi||2  =  0  and  thus  the  output  power  becomes 


E  [ \Vk\ 2]  =E  [|cfesu[/c]  +gifw[fc]|2] 
=  as2||cdes||2  +  a2. 


Sq  —  J(Bqo)  —  aqoGq  HdesHdesGgago  +  Cw- 


The  best  coefficient  a;;o  for  0  < 
sented  by  a  vector  with  q  entries: 


i  <  q  can  be  repre- 


Then  Si  =  a^G^HdesH£sGgagi  +  cr2.  Since  Sq  > 
a^GgHdesH£sGgag  +  a2  for  any  aq,  Sq  >  Si .  ■ 

As  will  be  discussed  in  Section  IV,  the  number  q  of  combined 
eigenvectors,  or  equivalently  the  dimension  of  the  signal  space, 
L  +  M  —  P,  should  be  estimated  by  using  an  information  the¬ 
oretic  criterion  such  as  the  minimum  description  length  (MDL) 
[14],  [15], 


C.  Choice  of  Decision  Delay 

For  simplicity,  we  have  so  far  considered  the  case  of  zero 
delay  6  =  0.  In  practice,  however,  we  have  to  set  the  decision 
delay  appropriately  because  the  delay  affects  the  performance. 
We  expect  that,  as  in  blind  equalization[16],  an  appropriately 
chosen  delay  will  improve  performance. 

Let  us  first  consider  the  influence  of  the  delay  S  on  the  ef¬ 
fective  channel  impulse  response.  Depending  on  6,  the  effective 
channel  impulse  response  c  can  be  classified  into  four  cases. 
Case  1:0<6<L  +  M  —  P 
In  this  case,  Hibi  takes  the  following  form: 


Hibi  =  [H(:,  0  :  <5  —  1)  H(:,  S  +  P  :  L  +  M  -  1)] 

where  H(:,  i  :  j)  is  the  MATLAB  notation  denoting  the 
2th  to  yth  columns  of  H.  Then  the  interference  component 
becomes 


Cibi  =  [Co  •  •  '  C«- 1  Cg+P  '  '  '  CL+M-lY' ■ 

Since  H.i,,  has  still  full  column  rank,  Cibi  can  be  reduced  to 
zero  by  minimizing  the  MERRY  cost  function  Jg( g).  Con¬ 
sequently,  the  effective  channel  impulse  response  becomes 

c  =  [0  •  •  •  0  eg  ■■■  cg+P- 1  0  •  •  •  Of. 
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TABLE  I 

Influence  of  6  on  d6  and  l8 


5 

d.g 

h 

Case  1 

0  <8  <L  +  M - P 

S 

p 

Case  2 

L+M -P+l <8 <L+ M ~ 1 

8 

L  +  M  -  <5 

Case  3 

L  +  M  <5<N 

— 

0 

Case  4 

N  +  1  <  8  <  Q  -  1 

0 

8  -  N 

Case  2:  L  +  M  —  P  +  1<S<L  +  M—  1 

Hibi  =  H(:,  0:5—1) 

Cibi  =  [co  •  •  •  c$- i]r 

c  =  [0  •  •  •  0  eg  •  cl+M-i]T. 

Case  3:  L  +  M  <  8  <  N 

Hibi  —  H,  Cibi  —  c,  C  —  0. 

Case  4:  N+1<8<Q  —  1 

Hibi  =  H(:,8  —  N  :  L  +  M  —  1) 

Cibi  =  [cg-iv  •  •  •  cl+m- i]T 
c  =  [c0  •••  CS-N- 1  0---  0]T. 

Table  I  summarizes  the  influence  of  8  on  the  position  d$  that 
the  nonzero  part  of  the  effective  channel  impulse  response  starts 
from  and  the  length  Zg  of  the  nonzero  part.  The  length  lg  should 
be  long  as  possible  because  the  power  of  the  desired  compo¬ 
nent  at  the  equalizer  output  will  increase  as  this  length  becomes 
longer.  Since  we  can  see  in  Table  I  that  the  length  lg  is  affected 
by  the  delay  8 ,  the  choice  of  the  delay  is  important. 

Exhaustive  determination  of  the  delay  that  provides  the 
highest  output  SINR  would  impose  a  heavy  computational 
burden,  so  delay  is  better  determined  by  a  computationally 
simple  technique  like  that  used  in  [12],  Here  we  propose  a 
simple  technique  to  determine  the  delay  without  estimating  the 
channel  impulse  response  directly.  Our  heuristic  idea  is  to  base 
the  search  for  a  delay  on  the  observation  of  the  original  channel, 
instead  of  the  effective  channel.  Consider  the  autocorrelation 
of  the  received  signal  defined  by 

Pd  =  E[  rf_1+drQ_i+d],  (23) 

Taking  into  account  the  cyclostationarity  of  the  transmitted 
signal  Uk,  we  can  rewrite  this  function  as 

P-l+d  0 

E  IN  ||“,  0<d<M-P+l 

i=d 

M  9 

as  E  ||h[i]||  ,  M-P  +  2<d<M 
0,  *  M  +  1  <  d  <  N 

d-N-l  „ 

E  ||hH||  ,  N+l<d<Q-l. 

7  —  0 


The  quantity  pd  measures  how  the  channel  impulse  response 
concentrates.  A  delay  is  determined  as  follows: 

8  =  arg  max  Pd ■  (25) 

0<d<Q-l f 

The  best  delay  is  thus  determined  by  computing  the  energy  of  a 
part  of  the  channel  impulse  response  and  searching  for  the  point 
from  which  the  maximum  energy  part  starts.  Unlike  the  method 
used  in  [12],  this  method  is  robust  to  the  channel  noise  because 
the  noise  component  vanishes  in  pd„  The  validity  of  this  method 
is  shown  by  computer  simulation  in  Section  V. 

IV.  Implementation 

The  blind-channel  shortening  method  presented  in  the 
previous  section  can  force  the  IBI  component  at  the  equalizer 
output  to  be  zero  even  in  the  presence  of  noise.  Since  this  means 
that  the  channel  can  be  shortened  within  the  CP  length,  the 
method  can  perform  perfect  channel  shortening.  The  proposed 
method  can  be  summarized  as  follows:  After  the  matrix  Gq 
specified  in  (19)  is  formed  by  linearly  independent  vectors 
that  minimize  Jg(g),  the  vector  aq  that  maximizes  J(aq)  is 
computed.  Then  the  equalizer  parameter  vector  g  is  computed 
according  to  ( 1 8).  In  this  section,  we  present  batch  and  adaptive 
algorithms  implementing  the  proposed  method.  Before  the 
equalizer  parameter  vector  is  computed  in  the  proposed  algo¬ 
rithms,  the  decision  delay  is  determined  by  using  the  algorithm 
in  (25)  with  the  time-average  of  pd  over  Dq  blocks 

Bg  —  l 

*  J-  V 

Pd  73  /  ,  kQ-\-P—  l-\-d^kQ-\-Q— 1+d*  (“6) 

0  k= 0 

A.  Batch  Algorithm 

The  batch  algorithm  estimates  the  correlation  matrices  R* 
and  Rr  by  time-averaging  and  then  eigendecomposes  them. 

Step  1:  Obtain  the  time  average  estimate  of  R*  over  13  \  blocks 
as 

1  k=i 

where 

f[fc]  =  r[(fc  -  1  )Q  +  P  -  1  +  S\  -  r[(fc  -  1  )Q  +  Q-1  +  8}. 

Step  2:  Eigendecompose  R*  and  obtain  eigenvalues 

An  >  A|  >  •  •  •  >  Aip_i  and  corresponding  eigenvectors  g,, 

i.  —  0.1,...,  l.  p  —  1 . 

Step  3:  Estimate  the  dimension  of  the  noise  subspace  of  R5  by 
using  the  MDL  signal  rank  test  [15] 

q  =  Lp  —  arg  min  MDL[i] 

0<i<Lp— 1 

and  form  G(/. 


Pd  = 


(24) 
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Step  4:  Obtain  the  time  average  estimate  of  over  B 2  blocks 

as 

Rr  =  -^Y^r[k]rH[k]. 

2  k= 1 

Step  5:  Eigendecompose  G^RrGg  and  obtain  the  eigenvector 
aq  corresponding  to  the  maximum  eigenvalue. 

Step  6:  Obtain  the  parameter  vector  as  g  =  Gqa q. 

B.  Adaptive  Algorithm 


tational  complexity.  Nevertheless,  it  is  highly  desirable  that  it 
need  not  estimate  the  noise  subspace  dimension.  There  would  be 
more  sophisticated  algorithms,  but  this  simple  algorithm  works 
as  shown  in  the  next  section  and  is  meaningful  for  the  adaptive 
implementation  of  our  SINR  maximization  approach. 

We  now  consider  the  convergence  property  of  the  proposed 
adaptive  algorithm.  The  proposed  adaptive  algorithm  consists 
of  two  parts,  i.e.,  the  adaptation  of  Lp  gs  and  that  of  a.  Both  of 
them  are  based  on  the  minimum  (maximum)  variance  method 
with  norm  constraint  as  in  the  MERRY  algorithm  [11].  The  con¬ 
vergence  analysis  of  the  MERRY  algorithm  has  been  shown  in 
[11]  and  can  be  applied  to  ours  directly.  Consequently  the  pro¬ 
posed  algorithm  can  be  shown  to  be  global  convergent. 


In  the  adaptive  algorithm,  vectors  that  minimize  Jg  are  ob¬ 
tained  by  using  the  adaptive  MERRY  algorithm  to  update  Lp 
vectors  with  different  initial  values  in  parallel.  A  stochastic  gra¬ 
dient  technique  is  used  to  obtain  a^l . 

Step  1:  Initialization: 

Determine  step  gains  //  and  /t„  appropriately.  Let  e,  be  the  ?'th 
coordinate  unit  vector.  Set  the  variables  as  follows: 

gt [0]  =  ei+i,  i  =  Lp  -  IjLp  -  2, . . .  ,0 
a[0]  =  [oo[0]  ai[0]  •  •  •  olp-i[0]]T  =  [1  0  •  •  •  Of  k  =  1. 


Step  2:  Estimation  of  eigenvectors: 

1  +  (>]  —  r  [(k  —  1)Q  +  Q  —  1  +  <5] 


f[fc]  =r[(fc  -  1  )Q  +  P  ■ 

£i[k]  =gf[k-l]r[k] 
n, i  =  Si[k  ~  1]  ~ 

84  J  \\gl[k-i]-^*mk}\\ 


i  =  0, 1, . . .  ,Lp  —  1. 


C.  Computational  Complexity 

In  the  batch  algorithm,  the  computaitonal  complexityies  of 
each  step  are  0(L2p2Bi)  for  Step  1,  0{L3p3)  for  Step  2, 
0(Lp )  for  Step  3,  0(L2p2B2)  for  Step  4,  0(L2p2q )  for  Step 
5,  and  O(Lpq)  for  Step  6.  The  major  computational  cost  is 
to  compute  eigendecomposition  of  the  correlation  matrices. 
Its  total  computational  complexity  is  proportional  to  L3p3. 
Although,  the  complexity  of  the  batch  algorithm  is  high,  our 
preliminary  simulation  results  showed  that  its  performance  is 
better  than  that  of  the  proposed  simpler  adaptive  algorithm 
when  an  amount  of  data  is  limited.  In  the  adaptive  algorithm, 
the  computational  complexities  of  Steps  2  and  3  are  0(L2p2) 
that  is  higher  than  that  of  the  adaptive  MERRY  algorithm. 
Although  the  computational  complexity  of  the  proposed  algo¬ 
rithms  is  not  low,  computational  complexity  should  not  be  the 
major  concern  since  microprocessors  are  getting  faster  all  the 
time.  Most  mobile  cost  lies  in  RF  and  analog  signal  processing 
[17],  [18]. 


V.  Simulation  Results 


Step  3:  Equalizer  output  power  maximization: 


A.  Simulation  Conditions 


G[fc]  =  [so[k]  gi [k\  ■  ■  •  gLP-i[fc]] 

VkQ  =a H[k  -  l]Gff[fc]r[(fc  -l)Q-l  +  8] 

a[k-l]  +  paytQGH[k]v[{k-l)Q-l+8] 

a[fc]  = - - - 

a[k-l\  +  paylQGH[k\r[{k-l)Q-l+8] 

g[fc]  =  G[fc]a[fc]. 


Step  4:  Set  k  =  k  +  1  and  repeat  Steps  2  through  3  until 
convergence. 

In  the  algorithm,  there  is  not  any  procedure  for  making  the 
estimated  eigenvectors  linearly  independent.  One  might  be  con¬ 
cerned  that  the  vectors  could  be  the  same.  Interestingly,  in  our 
simulation  of  static  channels,  the  vectors  differed  from  each 
other  by  setting  their  initial  values  different  from  each  other. 
Moreover,  in  the  algorithm,  all  possible  Lp  vectors  are  updated 
at  the  same  time  and  the  estimation  procedure  of  the  dimension 
of  noise  subspace  is  omitted.  Indeed,  there  are  useless  vectors 
among  Lp  vectors  since  the  dimension  of  the  noise  subspace 
of  R,s  is  less  than  Lp.  This  would  slightly  increase  the  compu- 


In  this  section,  we  show  computer  simulation  results  illus¬ 
trating  the  performance  of  the  proposed  algorithms.  The  infor¬ 
mation  bits  were  mapped  onto  QPSK  symbols.  A  raised-cosine 
pulse  p(t)  with  a  roll-off  factor  of  0. 1  was  used.  The  pulse  was 
truncated  to  by  6 Ts.  The  performance  measure  was  the  SINR 
defined  as  follows: 


SINR  = 


E  [\c*ud*[k]\2] 
E[\yk-c%ud*[k}\2] 

°IM2 


-IMI2)  +  ^ 


(27) 


where  cj  =  \c,i  ■  ■  ■  Cd+p]T ,  Ud[fc]  =  [uk-d  ■■■  Uk-d-p]T , 
and  d*  =  argmax^  1 1 1 1 2 .  Thus,  Cd*  represents  the  signal  com¬ 
ponent  that  is  the  largest  energy  consecutive  P  +  1  samples. 
Average  SINR  was  obtained  by  averaging  the  SINR  over  100 
different  channels.  The  received  SNR  was  defined  by 


SNR  = 


E  [| rk  -  wk\2 

E[\wk\2} 


M 


E  l|h[i]||s 


i=0 

P°l 


(28) 
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Decision  delay:8  Number  of  combined  eigenvectors:  q 


Fig.  2.  Influence  of  6  on  the  SINR. 


Fig.  4.  Influence  of  the  number  of  combined  vectors  on  the  SINR. 


Fig.  3.  Autocorrelation  pd. 


Fig.  5.  SINR  performance  of  various  batch  algorithms. 


B.  Batch  Algorithm 

We  set  N  =  512,  P  =  16,  L  =  45,  and  p  =  2.  We  consider 
nine -ray  multi-path  channels  whose  overall  response  is 

t  ~  +  &2 p(t  ~  5 Ts) 

+  hp(t  -  10TS)  +  b4p(t  -  15 Ts)  +  b5p(t  -  20 Ts) 

+  bfij)(t  -  25 Ts)  +  b7p(t  -  30 Ts)  +  bHp(t  -  35 Ts). 

Channel  coefficients  bi  were  zero-mean  complex  Gaussian 
random  variables,  and  the  length  of  the  original  channel  im¬ 
pulse  response  was  M  +  1  =  41.  Unless  otherwise  stated, 
SNR  =  50  dB,  B0  =  10,  and  Bi  =  B2  =  1000. 

We  first  examined  the  validity  of  pq  for  the  determination  of 
a  delay.  For  a  given  channel,  SINR  is  plotted  as  a  function  of  S 
in  Fig.  2,  where  we  can  see  that  the  SINR  strongly  depends  on 
6.  Fig.  3  shows  the  autocorrelation  pq.  As  can  be  seen  in  these 
figures,  the  SINR  tends  to  be  high  when  pq  is  high. 

We  then  examined  the  effect  of  combining  linearly  indepen¬ 
dent  vectors  in  the  noise  subspace  of  Rg.  The  average  SINR 
is  shown  in  Fig.  4  as  a  function  of  the  number  q  of  combined 
vectors.  In  this  case  the  dimension  of  the  noise  subspace  of  R* 
is  theoretically  21.  As  expected,  the  SINR  increases  monoton- 
ically  as  long  as  q  <  21.  This  result  shows  that  the  proposed 


algorithm  is  superior  to  the  MERRY  algorithm,  which  corre¬ 
sponds  to  q  =  1. 

We  next  compared  the  performance  of  the  proposed  algo¬ 
rithm  with  those  of  the  MERRY  algorithms.  We  had  various 
algorithms  use  the  decision  delay  determined  by  using  in  (25) 
with  the  value  of  pq  specified  by  (26).  Fig.  5  shows  the  average 
SINR  performance  of  various  algorithms.  Since  IBI  is  dominant 
in  high-SNR  regions,  the  performance  of  the  batch  MERRY  al¬ 
gorithm  with  p  =  1,  which  is  incapable  of  the  perfect  IBI  can¬ 
cellation,  is  very  poor.  Although  the  batch  MERRY  algorithm 
with  p  =  2  is  capable  of  canceling  IBI,  its  SINR  performance 
is  not  good  enough  high  because  of  noise  enhancement  espe- 
sially  in  low-SNR  region.  The  proposed  batch  algorithm  per¬ 
forms  best. 

An  original  channel  impulse  response  and  the  corresponding 
effective  channel  impulse  response  obtained  by  various  batch 
algorithms  are  shown  in  Fig.  6.  In  the  case  of  MERRY  with 
p  =  1,  residual  IBI  can  be  seen.  In  the  case  of  MERRY  with 
p  =  2,  although  IBI  could  be  canceled,  the  magnitude  of  the 
effective  channel  impulse  response  equivalently,  that  of  the  de¬ 
sired  signal  component  was  not  large  enough.  The  proposed  al¬ 
gorithm,  however,  could  cancel  IBI  completely  and  provied  a 
large  effective  channel  response. 

The  CP  length  could  be  insufficient  by  accident  or  design. 
SINRs  are  plotted  in  Fig.  7  as  a  function  of  CP  length.  From 


h(t )  =  b0p(ti)  +  bip 
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Original  channel 


Time 

Effective  channel  obtained  by  batch  MERRY  algorithm  with  />=t 


Time 


Effective  channel  obtained  by  batch  MERRY  algorithm  with  p= 2 


Effective  channel  obtained  by  proposed  batch  algorithm 


Time 

Fig.  6.  Example  of  magnitude  of  channel  impulse  responses  with  various  batch 
algorithms. 


5  10  15  20  25  30  35  40 

Cyclic  prefix  length:  P 


Fig.  7.  SINR  versus  CP  length. 


Fig.  7,  we  can  observe  that  the  performance  of  each  algorithm 
degrades  as  the  CP  length  decreases.  This  suggests  that  when 
CP  is  available  the  channel  shortening  is  much  preferable  to 
the  channel  equalization  that  corresponds  to  the  case  in  which 
the  CP  length  is  zero.  This  motivates  us  to  consider  channel 
shortening.  When  the  CP  length  is  short,  the  performances  of 
MERRY  with  p  =  1  and  the  receiver  without  an  equalizer  are 
poor  because  of  residual  IB1.  When  the  CP  length  is  long,  the 


Number  of  blocks:  B  ,  B 0 

Fig.  8.  Effect  of  the  number  of  blocks  on  SINR  performance. 

performance  of  MERRY  with  p  =  2  is  the  worst  because  of 
noise  enhancement.  The  performance  improvement  of  the  pro¬ 
posed  algorithm  over  the  other  algorithms  depends  on  the  CP 
length,  and  is  not  always  great.  When  the  CP  length  is  long,  the 
performance  of  the  proposed  algorithm  is  close  to  those  of  the 
other  algorithms  since  the  influence  of  IBI  is  insignificant  so 
that  the  performance  degradation  of  the  other  algorithms  is  not 
serious.  From  Fig.  7,  the  superiority  of  the  proposed  algorithm 
is  evident  when  the  CP  length  is  not  too  short  or  not  too  long. 

Next  we  consider  the  effect  of  the  number  of  blocks  on  per¬ 
formances.  The  results  are  shown  in  Fig.  8  where  B  \  is  the  same 
as  B2-  As  seen  in  the  figure,  the  proposed  algorithm  is  superior 
to  the  other  ones  and  require  more  than  100  blocks  for  satis¬ 
factory  results.  The  number  of  blocks  required  for  satisfactory 
results  would  depend  on  conditions. 

C.  Adaptive  Algorithm 

We  set  N  =  64,  P  =  8,  L  =  20,  p  =  2,  and  SNR  =  30  dB, 
and  we  considered  four-ray  multi-path  channels  resulting  in  an 
overall  response  of 

h(t)  =  b0p(t )  +  bip(t- +  b2p(t  -  5 Ts)+b3p(t  -  10TS). 

The  length  of  the  channel  impulse  response  was  M  +  1  =  16. 
The  step  gains  for  the  proposed  algorithm  were  set  to  p  =  2  X 
10-3  and  pa  =  10-2,  and  the  step  gains  of  MERRY  with  p  =  1 
and  MERRY  with  p  =  2  were  respectively  set  to  p  =  2  X  10-3 
and  p  =  2  X  10-3.  These  step  gains  were  determined  such  that 
the  IBI  component  in  the  effective  channel  impulse  response 
was  as  small  as  possible  and  SINR  is  as  large  as  possible  after 
2,000  iterations.  Average  SINR  for  a  given  channel  was  obtained 
by  averaging  the  SINR  for  10  different  ensembles. 

The  SINR  curves  are  shown  in  Fig.  9,  where  we  can  see  that 
the  proposed  adaptive  algorithm  can  achieve  a  higher  SINR  than 
the  adaptive  MERRY  algorithm.  Fig.  10  shows  an  example  of 
the  original  channel  impulse  response  and  the  corresponding  ef¬ 
fective  channel  impulse  responses  obtained  by  various  adaptive 
algorithms.  As  in  the  case  of  batch  algorithms,  the  proposed 
adaptive  algorithm  is  superior  to  the  conventional  algorithms. 
From  these  results,  we  could  show  the  possibility  of  adaptive 
implementation  of  our  SINR  maximization  approach. 
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Fig.  9.  SINR  curves  obtained  with  various  adaptive  algorithms. 
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Fig.  10.  Example  of  magnitude  of  channel  impulse  responses  obtained  with 
various  adaptive  algorithms. 


As  for  adaptive  algorithms,  the  application  to  dynamic  chan¬ 
nels  is  important.  At  this  moment,  the  proposed  adaptive  algo¬ 
rithm  is  not  aimed  to  cope  with  the  time  variation  of  channels. 
Hence,  the  algorithm  should  be  modified  to  improve  the  tracking 
performance  on  fast  time-varying  channels. 

The  SNR  in  both  simulations  is  rather  high  to  show  the  su¬ 
periority  of  the  proposed  algorithms  clearly.  We  confirmed  that 
the  proposed  algorithms  work  in  a  lower  SNR  as  well. 

VI.  Conclusion 

This  paper  proposed  blind-channel  shortening  algorithms 
based  on  both  batch  and  adaptive  processing.  Unlike  the 
conventional  algorithms,  the  proposed  methods  can  not  only 


remove  IBI  but  also  maximize  the  output  SINR.  We  have 
considered  the  impact  of  the  decision  delay  on  the  performance 
and  proposed  a  simple  method  for  determining  the  appropriate 
delay.  Our  simulation  results  demonstrate  the  superiority  of  the 
proposed  methodss  over  conventional  ones. 
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